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Summary 

Lift dominated pointed aircraft configurations are considered in the transonic range. 
These are treated as lifting wings of zero thickness with aspect ratio of order one. An inner 
expansion which starts as Jones’ theory is matched to a nonlinear outer transonic theory 
as in Barnwell’s earlier work. New expressions for the wave drag due to the equivalent 
body are derived. Some examples of numerical calculations for different configurations are 
presented. 

1. Introduction 

In 1946, R. T. Jones (ref. 1) published a paper giving a formula for the lift and 
induced drag of ’’low aspect ratio pointed wings below and above the speed of sound”. 
The work presented here, and earlier in the references cited below, represents an extension 
of Jones’ ideas to the transonic range. It is reassuring that, under suitable circumstances, 
Jones’ formula for the lift and induced drag not only continues to hold but is even valid 
for wings whose aspect ratio is order one. Under these circumstances, also shock waves 
and wave drag generally appear. 

The basic ideas of how this type of flow behaves are set out in the report of Barnwell 
(ref. 2). The principal result is that the lift produces a flow that looks, in the outer region, 
like the flow past an equivalent axisymmetric body. This physical effect shows up in the 
inner and outer expansions used by Barnwell. Cheng and Hafez used similar ideas to define 
the apparent body and general equivalence rule in a series of papers (ref. 3, 4). Cramer 
(ref. 5) also studied the problem (with zero thickness as is done here) and essentially 
verified the results of Cheng and Hafez. 

In this paper and ref. 6, we have also considered wing-like configurations with zero 
thickness and aspect ratio 0(1) as in fig. 1. Thickness effects can be incorporated relatively 
easily. Inner and outer expansions are defined in essentially the same manner as Barnwell 
although the asymptotic matching is carried out in a different way using an intermediate 
limit. Then wave drag associated with the outer expansion is considered. Several compu- 
tations and an optimization are carried out to show the effect of planform and longitudinal 
distribution of lift on the transonic wave drag. 


* This work was partially supported by the Air Force Office of Scientific Research under 
grant AFOSR 88-0037, and Rockwell North American Aircraft. 
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2. Basic Equations and Boundary Conditions 

The problem is studied in the framework of inviscid aerodynamics. Since entropy 
increases across the shock waves are of third order in the flow perturbation, the full potential 
equation can be used as a starting point. The flow is thus, to this approximation, isentropic 


- = (-)’ 

Poo \PooJ 

The potential equation is an expression of the continuity equation 

V • pq = 0, p = density, q = velocity 
$ is the velocity potential such that q = V$, 


( 2 . 1 ) 


( 2 . 2 ) 


(a 2 -^)$„ + a 2 (*„ + *„) = 2* I (* ! ,«„ + $,*„)+24. i ,$ I $ ! „ + ^4. ! , 1 , + ^4. SI (2.3) 
and a is the local speed of sound \J ~ip! P- The total enthalpy integral can be written 

(2.4) 


a 


U 2 M 2 2 k U 2 ’ 


U = free stream speed, = Mach number at infinity. 

The boundary condition of flow tangent to the surface can be written 


V$ • VB = 0 (2.5) 

on B(x, y , z) = 0, which defines the surface. 

We consider here an untwisted wing of zero thickness specified by an angle of attack a 
and a camber function m(x). The chord of the wing c = 1 and the span 2 b is 0(1). Thus 

B(x,y,z) = 0 = y- af(x) + 0(a 3 ) ( 0 < x < 1 | (2.6) 

\ z le < z < z LE ) 

where f(x) = m(x) — x and m(l) = 1. For a straight trailing edge at y = 0, the trailing vor- 
tex sheet lies in the plane y = 0, x > 1. The planform is specified by (±zle(x)),zle( 1) = b. 

Another boundary condition that must be satisfied is the “Kutta condition” at a 
trailing edge where the flow is locally subsonic. This condition implies that the pressure 
loading at a trailing edge is zero. In approximations, such as the inner expansion which 
follows, each term satisfies this condition. Another interpretation of this condition is that 
unphysical pressure jumps are not allowed in the inner solutions. 

3. Inner Expansion and Far Field (r* — > oo) 

The approximation in general is based on a — + 0, — + 1 and in the usual transonic 

way ( K — = transonic similarity parameter, fixed). ei(a) = parameter of outer 

expansion, defined later. In the inner expansion, the observer remains a fixed distance 
0(1) from the wing, and distances are measured from the wing surface. 
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The inner limit process thus has 


a — ► 0 (x, y *, z\ A'fixed) 


where 


y* — y — af(x) 0 < x < 1 
= y x > 1 

The form of the inner expansion for the potential is thus 


(3.1) 


z\a\ M«.) = U{x +atpi(x,y*,z) + a 2 ^2(x,y‘,2) + 0(<* 3 )} (3-2) 

The presence of the overbar denotes the possibility of logarithmic switchback terms 
introduced into the inner expansion for purposes of matching with the outer expansion. 
Anticipating the result, we note here 


2 1 , . . . 1 


¥>2(3, y*, z) = log — =</>22(x) + log -=y>2i(*) + <p 2 (x, y*, z) 
V e i v e i 

Note also that the velocity components of the inner expansion are 


(3-3) 


jj- = 1 + a< ?i x + a2 (V2x ~ /Vy ) + 0(a 3 ) 

77 = a<pi , + a 2 ip 2 * + 0(o: 3 )etc. 

U y y 

Substituting the assumed expansion into the full potential equation we obtain the equa- 
tions for the first two approximations (in divergence forms, as follows from (2.2)) and the 
corresponding conditions of tangent flow 


V*V=0, <pi .(x,0,z) = f'(x) 


0 < x < 1 
— ZLE < Z < Zle 


(3.4) 




<fz ,(x,0±,z) = f'(x)<p lx (x,0±,z) 

y 


(3-5) 


where 


V* = ( - — , — ) = inner transverse gradient 
v dy* dz ' 6 

V*ipi = inner transverse velocity perturbation 


V * 2 


^ (P_ 

dy* 2 dz 2 


= inner transverse Laplacian 
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The first equation (3.4) can be thought of as the Prandtl-Glauert equation of linearized 
theory but with M 0 0 close to one. The second equation (3.5) shows how, in the inner 
representation, either compression y>\ x < 0 or expansion <pi x > 0, provides an effective 
volume source and can cause stream tube divergence. We note the expression for the 
x-component of the mass flux vector. 

JJj = 1 - « 2 + i(VV,) 2 ) + 0(a 3 ) (3.6) 

The quadratic form in (3.6) is almost the RHS of (3.5); the additional term comes from 

V-( Vll VVi). 

The inner expansion is the driver of the entire procedure. But as is now shown it 
is not valid as r* = \/y* 2 + z 2 tends to infinity. On physical grounds, we would expect 
the transonic flow far away, which in general contains shock waves, to be described by 
an equation of mixed elliptic-hyperbolic type. The Laplace eqn. (3.4) is of course always 
elliptic. 

Further, we show that the second term a 2 y> 2 becomes much greater than the first 
as r* — ► oo. The general symmetry of the solution <^ 1,^2 is 

<Pi( x \ y* , z) = — <^>1 (rc; — y*,z ) odd Lifting, vortex sheet 

< P 2 ( x ‘, y *, z) = v? 2 (z; — y*z ) even Non-lifting, source 

We study now the behavior of these solutions near infinity in r* . The solution for 
ipi can be represented by a dipole sheet (or vortices) or, most directly, by the use of the 
complex variable 


f = z + iy* — r*e'^ (3.7) 

The complex potential for the wing, which is flat in a cross-plane x = const., and for the 
trailing vortex sheet is 


ft + ith = - tjp -zl E ( x )j 0 < x < 1 

= 1) {« - x > 1 

The transverse components of velocity perturbation V’Vi = (tq , ) are found from 


(3.8) 


ti>i 


0<I 


< 1 


= -im { 1 - 


V? 




(3.9) 


X > 1 


The first approximation to the pressure distribution and the lift can be found from 


Poo ~~ V<W ’ a lo 

— — 1 - a 7^i* H 




(1 - » 1 - 0(7 - 1 Vi, + 


(3.10) 
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The dominant term is a dipole potential where the dipole strength D\{x) is equal to the 
lift h(x) up to the station x. (cf 3.15) 


<pi(r*,d*;x) 


h(x)cos8* 1 
2t r r* 1 P ’ 


(3.20) 


An unyawed symmetric planform has been assumed. For more general planforms we can 
use 


ZLE 

Di(*)= J 

-*LE 


[v ? i]j/=o dz. 


The potential problem for <^2 can be thought of as describing the flow past a thin wing 
with thickness and a volume distribution of sources. Thus, the far field contains a source 
term like logr* and a particular solution due to the RHS. From (3.20) the RHS, has a 
term 


Thus 


7 + 1 


2,2 , '•* 2 
V'x +( ^i * + 


d* 


7+ 1 ,l'i \ 2 cos 2 9* 
2 { '2« ) r* 2 



&jp 2 1 1 7+ 1 d 1 + cos 20* 1 

dr* 2 r* dr* r* 2 dd* 2 4 dx \ 2ir j r* 2 r* 4 


(3.21) 


Taking account of the particular solution, the far field of <^2 is 


<P 2 (r*,0*-,x) = 


7 + 1 W 
4 (2tt) 2 


log 


2 r* + S 2 (x) log r* + g 2 (x) 


7 + 1 Ai'l 

8 (2tt)2 


cos 29 + 0( 



(3.22) 


Thus, there is a non-uniformity as r* — * 00 (since ~ -pr, a 2 <^2 ~ a 2 log 2 r*) 

roughly when r*log 2 r* ~ This shows the need for an outer expansion. An expression 
for the source strength S 2 (x) can be found from the boundary value problem for <p 2 , but 
g 2 (x) is undetermined from an inner problem. g 2 (x) must be found by matching with the 
outer nonlinear boundary value problem. The presence of shock waves in the outer flow is 
reflected in g 2 . 

4. Outer Expansion and Near Field (r — * 0). 

The first few terms of the outer expansion necessary to match with <pi, (p 2 are consid- 
ered. The limit process associated with this expansion is the typical transonic expansion 
necessary to give the small disturbance equation (ref. 7). The representative point runs 
to infinity as a — > 0, — ► 1. The limit process associated with this expansion has 

(x, y, z, K ) fixed as a — * 0 where (y, z) = \/e7(y, z ). 

The general form of the expansion for the potential is 
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$(x, y, Z] Moo , a) = U{x + ei(a)(f> x {x, y, z; K) + e 2 (a)<f> 2 (x, y, 5; 7v) + e 3 (a)<£ 3 (x, y, z; K) ■ • •} 

(4.1) 

where t\ , e 2 , e 3 are found from matching. In order to match, it is necessary to obtain an RHS 
term similar to that in (3.5), 2yi(<^ 2 ) r . <f > 2 can be made to match with the dominant dipole 
of inner (p x and this forcing term then appears in the RHS of the equation for </> 3 . <j) X is a 
switchback type of function necessary for matching and turns out to be the axisymmetric 
flow produced by an equivalent body of revolution. Some details are now shown. 


q x 1 <9$ 

¥ ~ U~dx ~ 1 + €l< ^ lx + e2< ^ 2x + €z ^ x + " ’ 

q y 15$ 3/2 , . 1/2 , . 1/2 , , 

Jj = U~dy = Cl ^y + e2e i 4>2 ~y + e h y + ' ' ‘ etc ' 

Jp = 1 + *i( K - (7 - l)<Ai x ) - ^ 2(7 ~ 1)<A2 X - e 3 (7 - l)^ 3x + 0(e 2 ) 
The full potential equation (2.3) takes the form 


{H-ei(RT -(7 — l)^ la .) — e 2 (7-l)^2 x -e 3(7 _1 )^3 a: • . . — l-2e 1 <f>i x -2e 2 <t> 2 x -2 e 3 <j> 3x . . .} 
{ti'fiixx + € 2^2 xx + € 3^3 X x +•••} + {1 + 0(ei)}{efV 2 <^>i -f- eie 2V 2 <^ 2 4- eie 3 V 2 <^> 3 } 

= 2{1 + • • ■}e?(^i.^i a; . + <^1 ~<£i x ~) + ' 

<42 ^ 

V 2 = 1 

dy 2 5z 2 

Thus, choosing eie 3 = e 2 we have 

(K - (7 + l)^i x )^i xx + V 2 <£i = 0 , 

(RT - (7 + l)<f>x x )<f>2 xx ~ (7 + 1)^2 X ^1 XX + V 2 </» 2 = 0 ( 4 . 4 ) 

0(eie 3 ) (K - (j + l)<j>i x )<j> 3xx - (*y + l)<f> 3x <f>i xx +'V 2 <f> 3 = (7 + l) 4 > 2 x <t> 2 xx (4.5) 

<j> 1 satisfies the usual nonlinear transonic small-disturbance (K-G) equation. <f> 2 satisfies 
its linear variational equation and <j) 3 a forced variational equation. All the equations are 
of conservation type and can be written 

(if* ir -h+i)^ )l+ v.(V*) = 0 (4.6) 

( AV, 2j — (7 + l)'f'i r '? 2 J ix + V • (V<^ 2 ) = 0 (4.7) 


(4.2) 

(4.3) 
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(Khx - (7 + 1)*,*,)* + V • (Vfc) = { -l±ll(tf x ) x (4.8) 

The near field behavior r —* 0 of these solutions, obtained from V 2 4> = RHS, is given by 

<f>\{x\f) — Si(x) log r + Gi(x) + 0 (r 2 log 2 r) (4.9) 


. , ~ /i\ D 2 (x)cos0 , ..(S'lDi) .,1.. 2 - l-i - 

h(xi r,0) = — — + (7 + 1)— 1 — cos ^{-r log 2 r- -rlogr + -} 

Zn r Zn 4 4 4 

+ ((7 + 1) -° - ^ - - K^l) cos 0 { r - log f - j) + £>(r 3 ) 


(4.10) 


4>z(x\,r,9) + 


+ 


^ — 4 ~ j 2 2 ^ - 5 * * * * lo S 2 f + S 3 (x) log f + G 3 (x) - 

(7 + l ) 2 (S'^DJDJ )' bg3 . + 0( . 2 [og2 ?) 


(7 + 1 ) D)D)’ 

8 (2ir) 2 


cos 20 


16 


(2tt) 2 


(4.11) 


where r = y'' j / 2 + z 2 , 0 = tan -1 | = tan -1 | 

The source strength S 3 and doublet strength D 2 are found by matching with the inner 
solution. The source strength Si is found in a special way in the matching. The functions 
Gi(x), G 3 (x) are found when the boundary value problems defined by the singular behavior 
as r — ► 0 in (4.9, 10, 11) are solved (numerically). 

5. Asymptotic Matching 

A matching limit, intermediate to the inner and outer limits, is defined by a class of 
functions 77 (a) such that << 77 (a) << 1. A coordinate 

r, = 7j(a)r (5.1) 


is held fixed in this limit. Thus 

r v 

r = — — *• 00 , 

V 

In the intermediate limit the representive physical radius again runs to infinity as a — ► 0, 
M 0 o — 1 1 but not as fast as in the outer limit. For matching, the inner and outer limit 
expansions must read the same in the intermediate coordinate. Thus 



ac^i + o 2 (f ? 2 + •••<-> < j >\ + 62^2 + € 3 </> 3 + • • • 

(«-> denotes “matches to”) 

Note that in the matching 
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(5.2a) 


and 


, _ , \/c7r, 

log r = log 


1 


= log -2- - log -j= 

*1 \Ai 


log 2 r = log 2 — - 2 log — log — \= + log 2 —j= 
V V yci V e i 


(5.26) 


Note also that 


cos 6* 

r* 


V ~ <*/( x) 
y 2 + z 2 - 2ayf(x) 


cos 6 

r 


+ <xf(x) 


cos 2 9 
r 2 


+ ... 


(5.3) 


Writing these out using the near field expansions of this section and the far field expansion 
of the previous section we have 


a( 


/l ^ S - + ■•■) + a 2 log 2 -^=v?22(*) + a 2 log -j=z<p2i(x) + a 2 { ^4 ^ (2V) 2 l<5g2 ~n 

+ S 2 (x) log — + 0 2 (z) - — t — cos + ^K 7 ? 2 )} + • 




|si(x)(log ^ - log -^=) + Gi(x) + 


{ D 2 (x) cos 9 
2xv^r,/, + "' 


+ • 


+‘3 { (W ^ - ■ 2 log i log . -L + log 2 ^)+ 


(2tt) 


_r 

8 




+S 3 (x)(log ^ - log -j=) + G 3 (x) - (T 8 1} °2 <=” 2« + • • ■ } (*•*) 

Comparison of these two expansions shows that they match in an intermediate region 
with the choices 


J2_ 


a 






D 2 M = ii(i) 

« . . _ (7 + 1 ) Dip; _ (7+1) w 

lV J 2 (2? r) 2 2 (2? r) 2 

S3(x) = 5 2 (x),G3(x) = 92 ^) 


(5.5) 


In summary 



e 2 = ol 2 log 1/2 —j=, 
V e i 


€3 = a 


2 


e x (a) is defined implicitly by the relationship above. S x (x) is chosen by an internal switch- 
back in the outer expansion. The switchback functions in the inner expansion are 
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(5.6) 


, x 1 „ , % (7+1) « 

= - jSl(*) = 

<p 2 i(x) = G^x) - S 3 (x) = Gj(x) - S 2 (x ) (5.7) 

The principal physical result of the matching is the source distribution for the apparent 
body that generates the first axisymmetric outer potential r) 


o, r , (7+1) W 

Sl(l) - 2 ( 2 ff ) 2 


(5.8) 


This body depends only on the longitudinal distribution of lift h(x) . A correction axi- 
symmetric flow is provided by the source S 3 

S 3 (x) = S 2 (x) (5.9) 

which generates the axisymmetric part of <f> 3 (x, r,0).<f>z(x, f, 0) can be decomposed into 


using the form of <j> 2 


Then for fl 3 (x, r) 


<£ 3 (x, r, 0) = ft 3 (x, r) + ip 3 (x, r ) cos 20 


<f> 2 (x, r, 0) = ip 2 (x, r) cos 0 


(5.10) 


(5.11) 


(K- 7 + l**)n»„ - (7 + m^ixx + ^ 


(5.12) 


By considering the omitted terms, an overlap domain can be shown to exist for match- 
ing to this order. Also consideration of higher order terms in both expansions shows that 
the matching can be continued. Thus, the outer expansion reads 


where 


$ = U{x + or 2 log — r, 0) + 0(a 4 )} 

\Ai 

<j> = <f>i(x,r) + - — 4> 2 (x, r, 0) + - — <j> 3 (x,r,0) 

log 7 lo s^7 


(5.13) 


V*i 

It can thus be noted that the collection of terms <j>i,4> 2 ,4>3 which can be computed 
individually satisfy together the small-disturbance (K-G) equation 


(A'-( 7 + l)^** + V 2 ^ = 0 


1 


log 3/2 -1 


(5.14) 


7^ 
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6. Wave Drag 

There is of course induced drag associated with the trailing vortex system; the drag 
in dominant order, associated with <p lt is just that of Jones’ theory. From the point of 
view of induced drag, the wing considered here, which is flat spanwise, is an optimum. The 
spanwise circulation distribution (cf 3.14) is elliptical. The wave drag is connected to the 
shock wave system in the outer flow field. It could be calculated from the entropy increase 
in the wave system. 

For small disturbances to a free stream we have the result for the wave drag D w 

D w = po 0 T 00 c 2 J j [S ] s dydz + --- (6.1) 

where [S] a = jump in specific entropy across a shock. The integral is taken over all 
the shocks in the system. Using the expression for the entropy jump in transonic small 
disturbance theory (cf. ref. 7, cf. p 165 ff. for a discussion of wave drag) 

D w = - Poo U 2 log 2 4=c 2 f°° rdf P" [4> x ] \d9 + • • • (6.2) 

1Z Jo Jo 

Consider the differential conservation form associated with (5.14) 

f (k& - - |(V 2 </>) 2 ) + (fhtxh + \(4>*4>e)e = 0 (6.3) 

Integrating this divergence form over all space outside a small cylinder 


(— oo < x < oo) f c — ► 0 

around the x-axis enables the entropy jump of (6.2) to be related to radial momentum 
flow. (6.3) is not conserved across shocks so that shock jumps appear such as 





Let D w — • Then 


/ oo r/7 r roo r 

dx / d0<}> x 4>r\r c — .hm / rdr / 

-oo Jo xl — ttJr<! Jo 

If we consider the dominant term in (6.4) 

D Wl = dx d9 fa ) 

Jo Jo r ^° r 

From (4.9), <f>\ x = S' x (x) log r + Gi(x), = ^4^ and usin S 

jf 


de 


(W 


\X2 


(6.4) 


(6.5) 
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in which Si(l) = Si(0) = 0, we have 

D Wl =—27 r f Si(x)G , 1 (x)dx = 2v f S' 1 (x)G 1 (x)dx ( 6 . 6 ) 

Jo Jo 

This drag formula is exactly that of a slender body in transonic flow (cf ref. 7, pl61). 
Higher order terms in the drag formula can be found. 

7. Applications and Remarks 

Several applications have been made of the theory in its present form. For flat wings 
/' = — 1 . /j (x) is given by (3.15). The effective source strength for the equivalent body 

Si(x) = ( c ^ e( l n - 6-5). Equation (4.3) 

(K -( 7 + l)<f>i x )4>i xx + h~~ + -~<t> i~ = 0 ( 7 -l) 

is solved numerically with a small-disturbance code using (4.9) as the boundary condition 
for various K . The dominant term of the wave drag coefficient Cdw is calculated from eqn 
(6.6). The results are plotted as Cdw vs for two different angles of attack in fig. 2. 
Substantial drag due to lift is evident. The planform shape and the distribution of h(x) 
which is typical appears in fig. 3. 

Another set of calculations incorporates a parabolic body of revolution (thickness 
ratio .057) and adds the source strength of this body to Si(x). A series of planforms with 
semi-span zle(x) given by 

*«(.) = (7-2) 

and shown in fig. 4 was considered for various pi, = .995, a = .2 rad. The idea is to 
optimize the L/D figure of merit C Dw/(AR = aspect ratio) by a choice of planform. A 
minimum drag occurs for pi = 2.5. The planform shape and curve of C Dw vs. [i appears in 
fig. 5. Also shown in figs. ( 6 a - 6 c) for pL = 1.2,2,10 are isobars which make evident the 
shock wave which occurs. The wave drag for small /i is large because of the small sweep 
and for large pi because of rapid changes of l\ near the wing tip. These preliminary studies 
are meant to show the relative effectiveness of various planforms. 

It would be very useful to extend this work to give efficient ways of calculating the 
higher order terms in the wave drag. It is also possible to incorporate the effects of wing 
thickness S ~ (a 2 log -^-) into the formulation in a more systematic way. A first step in 

this direction is given in ref. 6 . 
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POINTED TRANSONIC WING 



Figure 1 Pointed Transonic Wing 


DRAG RISE DUE TO LIFT CHARACTERISTICS OF 
MODEL FIGHTER PLANFORM 

100x50 GRID 
SINE TIP FAIRING 



Figure 2 Drag Rise Due to Lift Characteristics of Model Fighter Planform 
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LIFT LOADING OF MODEL WING WITH 
SINE TIP FAIRING 



“a" FAMILY OF WING BODIES IN WHICH SEMISPAN 

UX-xM 

s = r — 

M-1 
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DRAG/ASPECT RATIO FOR p W+B FAMILY 
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